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Abstract 

We describe a light-cone wave function formalism for hadropro duct ion of heavy-flavors via 
gluon-gluon fusion. The developed approach is well suited to address coherence effects in 
heavy-quark production on nuclei at small X2 < 0.1- A~^^^. Nuclear attenuation in hadropro- 
duction of heavy-flavors is found to be similar to shadowing effects for heavy-flavor structure 
functions in deep-inelastic scattering from nuclei. Nevertheless, remaining differences in the 
theoretical formulation of both processes imply corrections to factorization theorems. 



1 Introduction 



Heavy-flavor production in hadronic interactions at high energy is one of the traditional 
apphcations of perturbative Quantum Chromodynamics (pQCD). Although perturbative 
calculations are available now in next-to- leading order pQCD, and seem to be in reasonable 
agreement with experimental data ([1] and references therein), the production of heavy 
quark-antiquark (QQ) pairs has certainly not exceeded its potential in reveahng precious 
details of strong interaction physics. One of the missing ingredients is the quantitative 
understanding of nuclear attenuation effects, which offer a possibility to learn about the 
propagation of colored QQ pairs through the nuclear environment. This matter is strongly 
connected to effects of "color-transparency" and "color-opacity" , which are heavily discussed 
issues within a large variety of processes, ranging from quasielastic electron scattering to 
heavy quarkonium production (see e.g. [2,3] and references therein). 

In leading order pQCD the cross section for heavy quark-antiquark production via the 
collision of hadrons hi and /i2 at a squared center of mass energy s, can be written as 

Here fl^{xi,iJ,^) and fl^{x2,iJ>^) are the densities of partons "i" and "j", carrying fractions 
xi and X2 of the light-cone momenta of the colliding projectile and target respectively. The 
partonic subprocess i + j — > Q + Q is described by the cross section aij{xiX2S, n^) and 
requires a squared center of mass energy X1X2S > 4:mQ, where mq is the mass of the heavy- 
quark. Factorization in parton densities and hard partonic subprocesses is carried out at a 
typical scale //^ ~ 4mQ. 

In this work we will be concerned with nuclear effects in inclusive hadroproduction of 
QQ pairs which carry a high energy in the laboratory frame, where the target nucleus is 
at rest. In particular we consider processes at small target light-cone momentum fractions 
X2 0.1 • A"^^^. Furthermore we will restrict ourselves to moderate xp, which is deflned as 
the fraction of the projectile momentum carried by the QQ pair as seen in the center of mass 
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frame. Thus we concentrate on the kinematic domain where QQ production is dominated 
by the gluon fusion subprocess g + g ^ Q + Q and neither annihilation of hght-quarks nor 
excitations of higher-twist intrinsic heavy-quark components of the projectile wave function 
[4] are of importance. The pQCD cross section ([l|) is then proportional to the density of 
gluons in the beam and target. In the case of open-charm production at typical Tevatron 
energies (s = 1600 GeV^^) this implies that we focus on the region < xp < 0.5. 

We develop a description of heavy-flavor hadroproduction in the laboratory frame in 
terms of light-cone wave functions (LCWF) [5]. We find, that the underlying QCD mech- 
anism of QQ production can be viewed as diffractive dissociation of projectile gluons into 
QQ pairs. The resulting formalism greatly resembles the description of coherence effects in 
diffractive scattering on nuclei. Although there is much similarity between nuclear shadow- 
ing effects in hadroproduction of heavy-flavors and shadowing of heavy-flavor contributions 
in deep- inelastic scattering (DIS), they are nevertheless different and cannot be cast in 
the form of modified nuclear gluon distributions which are unique for all hard scattering 
processes. 

We apply our formalism to open-charm production and find only small nuclear atten- 
uation effects in the kinematic region considered. For the conventional parameterization, 
(Tec oc A", where A is the nuclear mass number. Our results correspond to an exponent 
a ~ 0.99. The corresponding experimental situation is not absolutely conclusive yet (for 
a comprehensive review on open-charm production see [4,6]). In early experiments [7,8] 
open-charm production was estimated indirectly from the observed yield of prompt leptons, 
rather than by a direct observation of produced charmed particles. The strong nuclear effects 
found in these experiments, a ~ |, led to the suggestion of large intrinsic charm components 
present in the respective hadron projectile (see [4] and references therein). More recent ex- 
periments directly observe Z)-mesons and obtain a ^ 1 , which is in good agreement with 
our results (^((xf) = 0.24) = 0.92±0.06 in the WA82 experiment with a 340 GeV tt" beam 
[9]; a(0 <xf < 0.5) = 1.00 ± 0.05 in the E769 experiment with a 250 GeV 7r± beam [10]). 

Our presentation is organized as follows: In Section 2 we analyze the space-time picture 
of heavy-quark production at small X2. In Section 3 we develop a LCWF description for 
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hadroproduction of heavy-flavors on free nucleons. The central result in this section is the 
representation of the inclusive QQ production cross section in terms of the cross section for 
the scattering of a three-particle, color singlet, Q-Q-^uon state from the nucleon target. 
We also comment on the cancellation of comover interactions. We extend our discussion to 
nuclear targets in Section 4 and apply our formalism explicitly to open-charm. Finally we 
summarize and conclude in Section 5 and indicate the difference between nuclear attenuation 
in inclusive open-charm and exclusive charmonium production. 

2 The space-time picture of QQ production 
at small X2 

In the laboratory frame the quantum-mechanical mechanism of nuclear attenuation in 
hadroproduction of heavy QQ pairs at small X2 and moderate is similar to the mecha- 
nism of nuclear shadowing in deep-inclastic lepton scattering at small values of the Bjorken 
variable x — Q"^ /2Mu (with Q^ — —v^ -|- the squared momentum transfer and M the nu- 
cleon mass). In DIS the absorption of the exchanged virtual photon on the target proceeds, 
at X < 0.1, via the interaction of quark-antiquark fluctuations present in the photon wave 
function. These qq pairs are formed at a typical distance Iq^ from the target. 

For a; < 0.1 the propagation length of the qq pairs exceeds the average micleon-nucleon 
distance in nuclei, Iqq > Rnn- Consequently, the quark pairs will start to interact coherently 
with several nucleons inside the target. This leads to shadowing (see e.g. [11,12]). If x 
decreases further 

the propagation length of the quark pairs becomes larger than the size of the nuclear target 
itself 

\q > Ra- (4) 
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Hence, the quark- ant iquark pairs will interact coherently with the whole target nucleus while 
their transverse size |r| is frozen [11]. Shadowing will be fully developed then. 

Similarly, at small values of X2 and moderate xp hadroproduction of heavy-quarks can 
also be treated within the laboratory frame in terms of QQ fluctuations of gluons, interacting 
with the hadron projectile. Their typical propagation length is 

where = xiEiab denotes the laboratory energy of the respective parent gluon, while Eiab 
stands for the beam energy. Evidently, if the propagation length of the QQ pairs is smaller 
than the average internucleon distance in the target nucleus, Iqq < Rnn, QQ production 
on nuclei will be equal to the incoherent sum over QQ production on individual nucleons 
within the target. In this regime the intranuclear interactions of the beam gluons reduce to 
rotations in color space and the flux of gluons is preserved Q. However, if the energy of the 
interacting gluons in the beam is large or equivalently X2 < 0.1 ■ A~^^^ is small, one obtains 
as before 

hQ>RA, (6) 

which means that the formation of QQ pairs takes place coherently on the whole nucleus - 
in close similarity to DIS at small x. 

Once produced, a QQ pair evolves into charmed hadrons after a typical formation (re- 
combination) length 
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I 

''b 

where m£, is the D-meson mass. The inequality between If and Iqq greatly simplifies our 
further consideration, as it implies that we may neglect the intrinsic evolution of a produced 
QQ pair during its propagation through the nucleus. 

At large xp ^ I the momentum of the projectile is transferred to a heavy QQ state 
collectively by several partons of the hadron projectile [4]. Such a QQ pair can therefore 
not be assigned to a certain projectile gluon. Consequently, the excitation of heavy-quark 



^Effective nuclear attenuation is possible though via the tension of color strings formed by color exchanges 
(predominantly on the front surface of the target nucleus) [13]. 
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components of the projectile at ~ 1 is quite different from the QQ production mechanism 
at moderate xp- In this work we will focus on moderate xp only, where the above mechanism 
is not important. 

Prom a diagrammatic point of view an accurate analysis of coherence effects requires the 
calculation of higher order diagrams, in which the incoming projectile gluon interacts with 
more than one nucleon. In standard pQCD this becomes a nearly impossible task. As for DIS 
at X < 0.1 • A~^/^, an analysis of coherence effects is greatly simphfied if one uses a LCWF 
formalism. This allows one to make explicit use of the fact that at high-energies, or small 
X2.I the transverse size of partonic fluctuations is conserved due to Lorentz time dilatation. 
A very important difference between hadroproduction of heavy-flavors and DIS is, that in 
DIS one starts out with a color-singlet photon, whereas in the former case the incident gluon 
carries color charge. Therefore, one might have expected to encounter an infrared divergent 
QQ production cross section, unless the interactions with comoving spectator partons of 
the incident color-singlet hadron were fully included. We shall demonstrate however, that 
the LCWF formahsm of ref.[5] can be generahzed to hadroproduction of heavy-flavors in an 
infrared-stable manner. 

3 Hadroproduction of QQ pairs on free nucleons 

In this section we develop a formalism for hadroproduction of heavy-flavors in terms of 
LCWF valid at small X2 < 0.1 and moderate xp- Applied to hadron-hadron collisions our 
results can be shown to be equivalent to the conventional pQCD parton model. However 
our new formalism becomes truly valuable, when it is applied to coherence effects in hadron- 
nucleus interactions. 

First we will consider the production of heavy QQ pairs through the interaction of a 
gluon from the projectile with the nucleon target. To flnally obtain the total production 
cross section we have to convolute our result with the gluon distribution of the beam hadron. 

In the LCWF technique one must carefully distinguish between "bare" and "dressed 
(physical)" partons. In our speciflc process the incident gluon will evolve after the interaction 
with the target into a state which contains a QQ Fock component. If this component is 
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projected onto a final QQ state to calculate the production amplitude of heavy-quark pairs, 
one has to bear in mind that the physical gluon in the LCWF formalism has a non trivial 
Fock decomposition which also contains a heavy QQ state. 

To leading order in the QCD coupling constant, as, the incoming gluon G has a rather 
simple Fock decomposition. It contains bare gluons, two-gluon states, and quark-antiquark 
components. As discussed in Section 2, at small X2 < 0.1 the propagation length of the 
leading order Fock states exceeds the target size. Therefore, the ^'-matrix, describing the 
interaction of the beam gluon with the target, will not mix states containing a different 
number of partons. Furthermore and most important, S conserves the transverse separations 
and the longitudinal momenta of all partons present in a certain Fock state. Consequently, 
it is natural to work with light-cone wave functions in a "mixed" representation, given by 
these conserved quantities. 

In this representation the LCWF of the physical incident gluon can be written as (wher- 
ever it shall not lead to confusion, we suppress the virtuality of the incident gluon kf) 

\G{kl)) = ^l-nQQ-n^\g) + ^ ^g{z, f)\QQ; z, r) + ^ ^(010 ■ (8) 

z,r 5 

Here ^g(-2, t) is the projection of the physical gluon wave function onto a QQ state in the 
{z, r )-representation, where the light-cone variable z represents the fraction of the gluon 
momentum carried by the quark and r is the transverse separation of the Q~Q in the 
impact parameter space. Z^^^(OIC) represents the light quark-antiquark qq and gluon- 
gluon gg Fock components. The normalization of the bare gluon state in the presence of 
the QQ, qq and gg components is given by 

{g\G) = ^l-nQQ-n^, (9) 

where 

riQQ = J^'dz J d^f\^G{z,n\', (10) 

determines the weight of QQ states in the physical gluon. Similar, denotes the normal- 
ization of the light-flavor qq and gg components. The QQ wave function of the gluon 
can be obtained immediately from the QQ wave function of a photon ^-y* as derived in 
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ref. [11]. The only difference is a color factor and the substitution of the strong coupling 
constant for the electromagnetic one: 



'■em 

as 



{[z' + (1 - zf]e'Kl{er) + mji^o^(er)} . (11) 



(27r)2 

Kq i are the modified Bessel functions, r = |r |, and 

^ = z{l-z)kl + ml. (12) 

— * 

After the interaction with the target nucleon at an impact parameter 6, the incident 
gluon state is transformed into 

\G]h) ^ S{h)\G]h) . (13) 

Here S{h) is the scattering matrix in the impact parameter representation. The familiar 
eikonal form [14] of the S'-matrix for the scattering of two partonic systems a and h is 



Sih ) = exp 



(14) 



with the one gluon exchange potential (eikonal function) 



V(b) = ^ f ^^ e^v{ik-h). (15) 

are color SU(3) generators acting on the individual partons of a and b at transverse 
coordinates hi and bj respectively. The effective gluon mass fig is introduced as an infrared 
regulator. In heavy- flavor production rriQ ^ fig, which ensures that our final results will not 
depend on the exact choice of fig. Furthermore, as we will show shortly, fig can be absorbed 
into the definition of the target gluon density. To the lowest non trivial order, we have to 
keep terms up to second order in V(b) only 0. 

Unitarity of the S'-matrix at a fixed impact parameter b simplifies the calculation of the 

inclusive QQ production cross section. For every value of b it relates the probability Pqq, 
^The fact, that the color generators in S do not commute is not important to this order in as, because 
our final results will contain the symmetric piece of the product of two generators only. 
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to detect a QQ pair in the final state with the probabihties Pq and P^, to find a physical 
gluon or light-fiavor qq or gg states: 

P^ + Pqq + P^ = 1. (16) 

Let us first consider Pg-, which is defined through 

P(. = - {GfN'\S{h)\G'N){GfN'\S{Tb)\G'NY . (17) 

N and A^' stand for the initial target nucleon and the final nucleonic three-quark states. 
The indices "f and "i" characterize the color charges of the final and incident gluon, where 
we average over the latter. Using the explicit form of the S'-matrix and applying closure for 
the final nucleonic states N\ we may simplify Eq. (^) in order a| to 

Pg = - y^{GfGfN\S{h)\G'G'N) . (18) 

We explicitly used the fact, that { — (T")*} are the antiparticle color generators, and re- 
grouped the gluon wave functions in Eq. (|1^) accordingly. Carrying out the color sums 
leaves us with the elastic scattering amplitude for the interaction of a color-singlet two 
gluon state {GG)i with the target nucleon: 

PG={{GG),N\S{h)\{GG),N). (19) 

Making use of the Fock state decomposition (H) of the dressed gluon, we obtain 

Pg = {l-nQQ-n^){{gg),N\S{h)\{gg),N) 

+ 2 j dzd^r\mG{z,r)\''{{QQg)iN\S{h)\{QQg)iN) 

+ 2 / di\^>m\{ig),N\S{hm9)iN) • (20) 



To the considered order in pQCD the contributions of different flavors in (|T6D and (^0]) do 
not interfere. We therefore suppress qq and gg terms in the following. In the color-singlet 
[gg)i state in Eq. (|20|) both gluons enter at the same impact parameter. Consequently this 
state has a vanishing color-dipole moment and 

{{gg),N\Smgg)rN) = l. (21) 



From Eqs. ( pO| , |2TD and the explicit form (|1^) for n^g, we obtain 

Pqq = 2 [ dzd'r\^fGiz,r)\'{iQQ9)iN\ \l - 5(6)1 \iQQg)iN). (22) 



Note that Pqq takes the form of a profile function for the elastic scattering of a {QQg)i 
state on a nucleon. After integrating over the impact parameter, we finally obtain the cross 
section for the inclusive production of heavy-flavor quark-antiquark pairs in gluon-nucleon 
collisions: 



(t{GN^QQX) = 2 I dzdPr\^>G{z,r)\^ I d%{{QQg)iN\ I - S{h) \{QQg)iN) 



j dzd'^r\^G{z,r)\^a^{fi,r2) = (asAf) 



QQg^ 



(23) 



where 



2 d%{{QQg),N\ l-S{h) \m9)iN) 
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(T(ri) + (T(r2) - ^a{r) 



(24) 



As indicated a^N is the total interaction cross section for the scattering of a color-singlet 
three-parton state {QQg)i from a nucleon [5]. The g-Q and g-Q separation in the impact 
parameter space are labeled fi and r2 respectively. Their difference r = ri — r2 is identical 
to the Q-Q separation. The distances ri^2 are related to the light-cone momentum fractions 
carried by the Q and Q respectively: 



r2 = —zf, fi = [1 — z)f 



(25) 



The cross section a{r) for the interaction of a QQ color-dipole of size r = |r | with a nucleon 
reads in Born approximation [11] 

_ 16 r d^k 



-q;| (l — J-'2{k)^ (l — exp{ik ■ 



(26) 



3 J (A;2 + /.2)2 

where J^2{k) = (A^|exp(iA; • (pi — p2))|A^) is the two-quark formfactor of the nucleon. In 
general, the dipole cross section is related in the Leading-Log(Q^) approximation (LLQA) 
to the gluon distribution of the target [5,16,17]: 



cr(r) a(x2, r) 



TT 



r'^as{r) 



A. 



X2g{x2, k2 ^ ^) 



(27) 
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where ^ ~ 10 [18]. The exphcit X2 dependence in ( p7| ) resuhs from higher order Fock 
components of the incident physical gluon, e.g. QQg states, and is determined through the 
BFKL evolution equation (for a detailed BFKL phenomenology of the dipole cross section see 
[15]). The gluon distribution in ( pTl) absorbs the dependence on the infrared regulatization 
/ig in the same manner as in the conventional QCD analysis of hard scattering processes. 
The most important point for our further discussion is the color transparency property, i.e. 
the oc dependence of the dipole cross section in (p7|). 

When cast in the form (p3|) , the heavy- flavor production cross section resembles heavy- 
flavour contributions to real and virtual photoproduction. While the three particle cross 
section ( [2^ ) enters in the former, the dipole cross section a{r) is present in the latter. (See 
also Eq. ([11]) for the relationship between the corresponding wave functions.) Both processes 
have a similar dipole-size structure, and we briefly recapitulate the results of refs.[5,ll]: The 



QQ wave function (11) decreases exponentially for r > 1/e. We may therefore conclude that 
for kf < ArriQ small transverse sizes are relevant for the QQ production process: 

r^rlrl<^. (28) 

For highly virtual gluons, kf ^ 4mQ, the QQ production cross section ( p3D vanishes ^ 1/kf 
due to the color transparency property of the dipole cross section. (In view of the analogy to 
high energy photoproduction, this behavior can be understood as the counterpart of Bjorken 



scaling.) Combining equations (^^ ^), we therefore find 

a{GN QQX) oc X2g{x2, ~ 4mJ) . (29) 

The color transparency property ( pT] ) of the dipole cross section has important impli- 
cations. First note, that ^qq, the probability that a physical gluon is realized in a QQ 
Fock component, has a familiar logarithmic ultraviolet-divergence at r ^ 0. Nevertheless, 
in the LCWF formalism the heavy-quark production cross section a{GN QQX) is ul- 
traviolet finite, since the color transparency property of the dipole cross section ( P?] ) leads 
to a rapid convergence of the r-integration in the cross section (pB]). The physical reason 
for color transparency to play a role in this process is, that the small size QQ pairs cannot 
be resolved by interacting t-channel gluons with wavelengths A > r. They are therefore 
indistinguishable from bare gluons and cannot be excited into QQ states. 
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The second remarkable feature of the LCWF formahsm is that a{GN QQ X) is 
also infrared-stable. This is not the case for the total gluon-nucleon cross section, being 
(logarithmically) infrared-divergent for fig 0. Only the interaction of two color-singlet 
hadrons yields an infrared-finite total cross section. This fact is due to the cancellation of the 
incident gluon color charge through the color charge of comoving spectator partons belonging 
to the color-singlet hadron projectile. However, soft t-channel gluons with k ^ cannot 
resolve the QQ component of the physical beam gluon and therefore cannot contribute to 
heavy-flavor production. We consequently obtain an infrared-stable QQ production cross 
section p3| ) already for an isolated, colored, gluon projectile. 

A simple consideration shows, that the QQ production cross section (|23|) is not affected 
by the presence of spectator partons in the beam. Consider for example the case of one 
spectator "parton" s, which complements the incident gluon to a color-singlet state {Gs)i. 
The spectator parton is assumed to scatter elastically only, i.e. without dissociating into a 
QQ pair. Generalizing the above analysis, we finally end up with 



To obtain the inclusive production cross section we have to sum over the colors "/" and 
"d" of the final gluon and spectator states respectively. This yields a color-singlet two 
gluon state, as well as a color-singlet spectator ss state. Since both partons s and s of 
this singlet state enter at the same impact parameter, their color-dipole moment vanishes 
and they decouple from the interacting gluons. This leaves us with our previous result 
given in Eq. (pS]). A similar cancellation takes place for an arbitrary number of spectators, 
as long as we neglect interference effects between QQ pairs belonging to different gluons 
of the beam hadron. However such a restriction to elastic rescatterings of spectators is a 
standard ingredient of pQCD calculations in the framework of the conventional LLQA. In 
this respect the cancellation of interactions with spectator partons can be traced back to the 
LLQA factorization of the LCWF, containing the heavy-quark component. Consequently in 
LLQA, the production cross section of QQ pairs in the kinematic region under consideration 
is determined by the gluonic subprocess G + N ^ QQ + X, i.e. the underlying process is 
diffractive dissociation of gluons into QQ pairs. 




(30) 
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We are now in the position to write down the heavy-quark production cross section for 
hadron-nucleon colhsions aQqih, N). For this purpose we have to multiply a{GN QQ X) 
from Eq. ( p3D with the gluon density of the incoming hadron projectile and integrate over 
the virtuality kf of the projectile gluon: 

da,,(h.N) ^ f (31) 

dxi J kf ologkf 

As mentioned above, the leading contributions to a{GN QQ X) result from the region 
kf <4:'mQ, where we approximate a{GN — > QQ X]X2,k\) ^ cx{GN QQ X]X2,k\ = 0). 
We therefore obtain 



dcrQQ{h,N) 



gi{xi,p? = Ami) ^('^^ ^ <5<Q X- xa, A;^ = 0) . (32) 



dxi 

Notice, that we have recovered the dependence on the gluon density of the beam as in the 
parton model ansatz (|1|), whereas (|29| ) yields the gluon density of the target. This shows 
the equivalence of the LCWF formulation of heavy-flavour production and the conventional 
parton model approach. 

Let us explore the above result for charm production in nucleon-nucleon collisions. To be 
specific we choose a beam energy Eiab = 800 GeV, as used in the Fermilab E743 experiment 
[19]. We employ the dipole cross section of ref. [15], which has been shown to yield a good 
quantitative description of the closely related real photoproduction of open-charm, and 
the small x proton structure function, as measured at HERA. According to the analyses 
in refs. [18,20], one obtains for transverse sizes ~ V'^c ^ precocious asymtotic BFKL 
behavior 



a{x2,r) ^ aB{r) y—J , (33) 

valid at X2 < Xq. The Born cross section ^"^(r) from Eq. (^Bf) is evaluated for fig = 0.75 GeV, 
A = 0.4 is the intercept of the BFKL pomeron and Xq = 0.03. The charm quark mass is 
fixed at rric = 1.5 GeV. The gluon distribution of the nucleon projectile enters at moderate 
to large values of Xi. We may therefore use the parameterization of ref. [21]. Because the 
g + g ^ Q + Q cross section decreases rapidly with the invariant mass of the heavy-quark 
pair, one finds X1X2S ~ Aitiq. We may now calculate '^'^""Jf^^^ from Eq. (^). In Fig.l we 
compare our results with the data from the E743 proton-proton experiment [19]. In the 
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kinematic region, where our approach is well founded < xp < 0.5, we obtain resonable 
agreement with the experimental data. 



4 Hadroproduction of QQ pairs on nuclear targets 

Although the conventional pQCD result for heavy-flavor production and the LCWF formu- 
lation are equivalent, there is an important advantage in the latter. Namely in equation p^) 
for a{GN QQ X) the LCWF of the QQ state and its interaction cross section factorize. 
This feature is crucial for the now following generalization to nuclear targets. 

The derivation of Eq. (|23|) was based upon the observation that at high energies, where 
the coherence length of the QQ fluctuation of the incident gluon is larger than the target 
size Iqq > Rtarget, the trausvcrsc separations of the beam partons are frozen. This led to a 
diagonalization of the scattering S'-matrix in a mixed {z, r) -represent at ion. In high energy 
hadron-nucleus collisions with Iqq >jRa, this diagonalization of the S'-matrix can evidently 
be achieved too. Consequently we obtain the QQ production cross section for nuclear 
targets by substituting a^A for a^N in Eq. (pS]). The cross section a^A for the scattering 
of a three-parton, color singlet, {QQg)i state from a nucleus with mass number A is given 
in the frozen size approximation through the conventional Glauber formalism (we suppress 
the parton separations r,f 1^2)'- 



a3A = 2l£bll- 1 - ^ asNTib) d^b 1^1 - exp 



A' 

(34) 



CTsNTib) 



Here b is the impact parameter of the (QQg) -nucleus scattering process, which must not be 

confused with the impact parameter of the ((5(5fi')-nucleon interaction in Section 3. T{b) 

stands for the optical thickness of the nucleus 

_, r+00 _, 
T{b)= / dzuAib.z) , (35) 

with the nuclear density nA{b,z) normalized to / d^fnA{r) = A. In LLQA the inclusive 
cross section for QQ hadroproduction on nuclei can then be written as 



daQQ{h,A) 



gi{xi, jj^ = 4ml) ^ ^'^ ^2, kl = 0) , (36) 



dxi 
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where 



a{GA^QQX;x2,kl) = 

2 J dz(fr\^G{z,r)\^ j d^fejl -exp [-i a3^T(6 )] | = (a3^)QQg . (37) 



In the multiple scattering series ( |3^ , |37| ) nuclear coherence effects are controlled by the 
expansion parameter 

A^'\ (38) 

Expanding the exponential in (|3^j37| ) in powers of ta one can identify terms oc r^. They 
describe contributions to the total production cross section, resulting from the coherent 
interaction of the (QQg) state with u nucleons inside the target nucleus. In leading order, 
z/ = 1, we remain with the incoherent sum over the nucleon production cross sections: 

a{GA^QQX) = j (fhT{h) a{GN ^ QQ X) 

= Aa{GN ^QQX). (39) 

This yields the conventional impulse approximation component of the nuclear cross sec- 
tion, which is proportional to the nuclear mass number. Effects of coherent higher order 
interactions are usually discussed in terms of the nuclear transparency 

rp _ (^QQih,A) 

In the impulse approximation = 1. The driving contribution to nuclear attenuation 
results from the coherent interaction of the three-parton {QQg) state with two nucleons 
inside the target: 

Ta = 1-\ / d%T\h) . (41) 



4 {(t^n)(. 

At kl -C 4mQ shadowing effects can easily be estimated, using the fact that |\E'g(-25 ''")P oc 
exp(— 2mQr) and cr(r) ^ Cr^. We further approximate cr^{ri,r2) ~ |cr(r), and obtain 



(42) 



'QQg ^'"Q 

For Gaussian nuclear densities with mean square charge radii Rch = 1.1^^^^ fm we find 



/ 



d=6r^(6)^^. (43) 
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Altogether, we finally obtain the following estimate for nuclear attenuation effects: 

At pre-asymptotic energies, when the coherence length Iqq < Ra, the shadowing term 
in ( ^11) will be proportional to the squared nuclear charge form factor F^/t), where k = 
1/Iqq = Mx2 is the longitudinal momentum transfer in the G + N QQ + transition 
[22]. This form factor quantifies the onset of nuclear coherence effects. 

Applied to open-charm production at current Tevatron energies (s = 1600 Gel^^) we 
have to remember that X2 is bound from below through the kinematic constraint 

XiX2 > ^ ^ (0.5 - 1) • 10-2. (45) 
s 

In the accessible region of small X2 ~ 0.01 we find at r ~ l/m^ typically C ~ 2.5 [15], which 
leads to 

1 - Ta ~ 6 ■ 10-^ ■ . (46) 

Approximated by the ansatz T4 ^ A°'~^, we get 1 — a ~ 7 ■ 10"'^, in agreement with more 
recent experimental results [9,10]. However, this also shows that data with high accuracy 
are necessary to observe nuclear attenuation of open-charm production. 

Although nuclear shadowing effects turn out to be small, it is interesting to consider 
their consequences for factorization theorems. Remember the close similarity between nu- 
clear shadowing in heavy-flavor hadroproduction and heavy-flavor contributions to nuclear 
structure functions. It mainly involves an exchange of the three-parton cross section a^^ in 



(|41|) with the dipole cross section cr(r), which differs in the considered process by a factor 
~ 2 (see the discussion following Eq. (0)). Indeed, comparable ~ 1% nuclear shadowing for 
the charm structure function of lead was found in ref. [16]. Following the analysis of ref.[ll] 
we note, that the shadowing effect in Eq. (|4l|) does not vanish at high virtualities of the 
projectile gluon kf ^ 4mQ, and consequently is a leading twist effect. (The same is true for 
nuclear shadowing in DIS.) Our formalism correctly reproduces the proportionality of the 
QQ hadroproduction cross section to the gluon structure function of the target nucleon, in 
agreement with the conventional parton model formulation as shown in Section 3. A simi- 
lar proportionality to the gluon structure function of the target nucleon holds also for real 
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(and weakly virtual < 4m^) photoproduction of open-charm on nucleons. It is therefore 
tempting to assume such a factorization also in nuclear production processes. However the 
three-parton cross section and the dipole cross section, entering hadro- and photoproduction 
respectively, are not equal. Hence the A-dependence of both processes cannot be described 
in terms of modified nuclear gluon distributions, which would be unique for all hard pro- 
cesses. The fact, that nuclear shadowing defies factorization, despite being a leading twist 
effect, was already encountered earlier in connection with DIS and Drell-Yan production on 
nuclei [16], and also in nuclear production of vector mesons [23] and jets [24]. 

A detailed discussion of the subtleties of factorization breaking goes beyond the scope 
of the present paper. Here we only wish to note, that in DIS and heavy-flavor hadropro- 
duction non-factorizable effects result from the contributions of QQ Fock states present in 
the projectile photon or gluon respectively. Furthermore, shadowing of QQ contributions in 
DIS must be reinterpreted as a modification of the nuclear sea quark distributions [5] . This, 
however, is not possible for heavy-flavor hadroproduction. Shadowing effects, which can be 
ascribed to nuclear modifications of gluon structure functions, emerge in DIS first through 
contributions from higher order QQg Fock states of the photon. For the scattering from 
free nucleons these lead to the conventional QCD evolution of the target gluon distribution 
function [5] included in Eq. (^). However in the case of nuclear targets the situation be- 
comes different as we will outline briefly. In DIS at large momentum transfers a QQg Fock 
state, which contains a soft gluon, is characterized through a small transverse separation of 
the quark-antiquark pair Tqq ~ Ij ^Q'^ + 4mQ, but a large separation r ^ Rg = 1 / fig of the 
gluon from both quarks (see refs.[5,25]). Note that the latter is independent of the photon 
virtuality Q^ and rriQ. A QQg Fock state of the virtual photon can therefore be treated as 
an octet-octet color dipole of size r, because the QQ pair acts as a pointlike color charge by 
virtue of r ^ ^^qq- Furthermore the squared wave function of the QQg Fock state factorizes 
into the square of the QQ wave function |^g'(z, rgg)!^, and the distribution of soft gluons 
around the nearly pointlike QQ pair: 

\^,{z,rQQ,Zg,r)\^ = ^as{rQQ)\^,,{z,rQQ)\^^^ (47) 

Here Zg is the momentum fraction carried by the soft gluon and $(x) = x^[Kf{x) + 
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xKi{x)Ko{x) + ^x'^Kq{x)]. Nuclear attenuation of QQg states will be controlled by the 
color octet dipole cross section <7s{r) = |cr(r) [5]. On the other hand, the contributions of 
QQg components to shadowing can be related to diffractive dissociation of virtual photons 
into large mass states ^ Q^, which is successfully described in terms of the so-called 
triple-pomeron coupling A^jp. At X2<xa = 0.1 A^^^ one finds [5,22,26] 

1 - ^ I^AsiP log (^) ~ 1.5 ■ 10-^AV3 log . (48) 

We used A^jp ^ 0.16 GeV^^ as determined in the photoproduction experiment [27]. 

Similar considerations are applicable to the excitation of QQg states in hadroproduction 
of heavy-flavors. Postponing detailed numerical calculations, we may only cite the main 
ideas: The underlying process will be diffractive excitation of incident projectile gluons 
Gi into quark-antiquark-gluon states Gi QQg2- Consequently the generalization of the 
nucleon (p3|) and nuclear (|37D production cross section will contain the interaction cross 
section of a four-parton color singlet state QQgig2, where the transverse size rgg ~ l/rnq 
of the QQgi component is much smaller than the separation r ~ i?g of the soft gluon g2. 
Furthermore they will involve the factorized three-parton wave function from Eq. (0). Since 
the size of the three-parton QQg state in DIS, as well as the size of the four-parton QQgig2 
state in hadroproduction is determined by the soft gluon separation r ~ 1/Rg, nuclear 
attenuation of QQg states will be similar in both processes. Only X2 has to be substituted 
for X in Eq. 

At Rg > Tqq the triple pomeron coupling becomes approximately independent of Q^ 
and quark masses [5]. Therefore we may reinterpret the nuclear attenuation of QQg Fock 
states as nuclear modifications of target gluon distributions. These are then unique for all 
hard processes. However in DIS the non-factorizable contributions of light flavor qq pairs 
dominate nuclear shadowing down to x < 10^^. Higher order Fock states are important 
at much smaller x only [5]. Hence the observed nuclear shadowing in DIS is not an ap- 
propriate place to determine nuclear modifications of gluon distributions. In open-charm 
hadroproduction the triple-pomeron contribution ( ^8]) will be relevant already at X2 ~ 10~^. 
Nevertheless possible attenuation effects will not exceed ~ (5-10)% even for heavy nuclei. 
This corresponds to an exponent 0.98 < a < 1 - a deviation from unity which is still below 



the current experimental accuracy. Of course Eq. ( ^S]) gives a crude estimate of the X2 
dependence of nuclear attenuation only. A more detailed calculation using the pomeron 
structure function of ref. [25] will be presented elsewhere. 

Finally, the estimate ( ^HD shall be applicable to nuclear shadowing in Drell-Yan and 
exclusive vector meson production on nuclei. In the latter case qq contributions to nuclear 
shadowing decrease asymptotically oc l/(<5^ + fny), but are numerically very large. The 
influence of color transparency [28] demands an extremely large momentum transfer for 
a substantial decrease of quark-antiquark qq contributions and the onset of the universal 
attenuation of Eq. (^8]). 

5 Conclusions 

The purpose of this paper was the presentation of a light-cone wave function approach to 
hadroproduction of heavy-quarks on nucleons and nuclei valid at small X2 < 0.1 A^^^ and 
moderate Xi? 7^ 1. We found that in the laboratory frame the mechanism of heavy-flavor 
production in LLQA can be described as a diffractive dissociation of projectile gluons into 
QQ pairs. The central result is the representation of the inclusive QQ production cross 
section in terms of the cross section for the scattering of a three-parton, color singlet, Q-Q- 
gluon state from the target. Although the incident gluon carries color charge, the light-cone 
formalism leads to an infrared-stable heavy-flavor production cross section, which is not af- 
fected by interactions with comoving spectators. Nuclear attenuation effects were found to 
be small ~ A/ {mQRAf' ■ We applied our formalism to open-charm production and obtained 
good agreement with more recent experimental data. Comparing nuclear effects in photo- 
production and hadroproduction of open-charm we observed, that despite beeing a leading 
twist effect, leading shadowing effects defy factorization. However nuclear attenuation of 
higher order production processes can be cast in the form of nuclear modifications of gluon 
distributions. 

Weak nuclear attenuation in inclusive hadroproduction of open-charm does not imply 
weak nuclear effects for exclusive charmonium channels. It should however be possible to 
extend the presented technique to the latter. As far as nuclear effects are concerned, the 
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difference between both should be quite similar to the well established difference between 
inclusive photoproduction of open-charm and exclusive photoproduction of charmonium 
bound states. In the first case nuclear attenuation is controlled by the cc dipole cross 
section at small distances r ~ l/mc- This leads to a weak nuclear attenuation [5,16]. On 
the other hand in exclusive production processes the typically transverse scale is given by 
the size of the charmonium bound states ^ ^/nT'c, with the consequence of much 

stronger nuclear effects [28] . A detailed treatment of exclusive charmonium production will 
however be presented elsewhere. 

Finalizing this paper we learned that B. Z. Kopeliovich is working on a similar analysis 
of nuclear shadowing effects in Drell- Yan processes [29] . 
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Figure caption: 

Fig.l Differential cross section of open-charm production in nucleon-nucleon collisions at 
Eiab = 800 GeV. The experimental data are from ref.[14]. 
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